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Using a standard cooling method for SU(3) lattice gauge fields constant Abelian magnetic field 
configurations are extracted after dyon-antidyon constituents forming metastable Q = configura- 
tions have annihilated. These so-called Dirac sheets, standard and non-standard ones, corresponding 
to the two (7(1) subgroups of the SU(3) group, have been found to be stable if emerging from the 
confined phase, close to the deconfinement phase transition, with sufficiently nontrivial Polyakov 
loop values. On a finite lattice we find a nice agreement of the numerical observations with the 
analytic predictions concerning the stability of Dirac sheets depending on the value of the Polyakov 
loop. 
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I. INTRODUCTION 

In lattice gauge theories the cooling method is 
used to remove short distance fluctuations in order 
to search for (approximate) classical solutions of the 
Euclidean field equations [J-Q- We consider this 
technique as a device 1 (like smearing or filtering 
based on low-lying modes of the Dirac operator) that 
may help to identify topological excitations generi- 
cally present in the sample configurations representing 
the zero-temperature (or thermal) ensemble of gauge 
fields [il]. 

Cooling studies of nonzero-temperature SU(2) lat- 
tice fields [H have identified as topological excita- 
tions both calorons with nontrivial holonomy [9rfllj or 
dyon-antidyon pairs which finally annihilate. Some- 
times this annihilation process provides a constant 
Abelian magnetic field called Dirac sheet (DS) , which 
turns out to be either stable or unstable under fur- 
ther cooling [12]. The stability is strongly correlated 
with the spatial average value of the Polyakov l oop 
(the holonomy) in the given stage of cooling. In [13( 
an explanation for this observation was presented. 

Some time ago we started cooling studies of SU(3) 
gluodynamics, applying the Cabibbo-Marinari proce- 
dure in the cooling mode for the standard Wilson 
action [l4j]. On plateaus characterized by values of 
the action within the range 0.5 — 1.5 times the one- 
instanton action Si ns t the emerging topological ob- 
jects turned out to be either calorons or anticalorons 
(dissociated or not dissociated into their three respec- 
tive dyon or antidyon constituents) or one or two 
dyon-antidyon pairs. Sometimes (similar to the SU(2) 
case) the annihilation process of a dyon-antidyon pair 
leaves behind a constant Abelian magnetic field. In 



the SU(3) case the structure of such Dirac sheets is 
somewhat richer than in the SU(2) case. Below we 
will describe their analytic construction following a 
seminal paper by Gerard 't Hooft [l5[. We will ex- 
pand the concept of marginal stability [l6l - fl8j to the 
SU(3) case. We shall find agreement between the an- 
alytically worked-out preconditions - in terms of the 
holonomy - for stability of the Dirac sheets in a finite 
volume on one hand and the numerical observations 
for Monte Carlo generated - and subsequently cooled 
- lattice gauge fields. 



II. DIRAC SHEET SOLUTIONS 

In lattice gauge theories usually periodic boundary 
conditions are applied for the gauge fields (by default, 
if no special needs suggest something else). Thus, 
the DS configurations, that can be obtained by the 
cooling procedure, are periodic as well. The simplest 
way, however, to present analytic solutions with a con- 
stant color-magnetic field on a hypertorus uses twisted 
boundary conditions [15]. In this case most of the 
structure of the solutions is absorbed into twists (the 
gauge transformations that the gauge fields acquire 
over the periods on a hypertorus). They look rather 
complicated and are even non- Abelian while the gauge 
fields themselves are rather simple. To have periodic 
solutions we should make clear that the twists can be 
removed by appropriate gauge transformations. The 
necessary condition for this is the commutativity of 
twists in different directions. Below we will apply this 
condition to find those solutions that allow to be made 
periodic. 

Discussing the special selfdual solutions, 't Hooft 



2 



was considering the general SU(N) case. The gauge 
field (x) and the field strength F ilv {x) are strictly 
Abelian while the twists are non-Abelian. The gauge 
field A^(x) is proportional to the diagonal traceless 
matrix ui — 2it diag(i, I, — k, — k) with positive 
integers I and k such that I + k = N 

F ll u{x) = -u}{a llv -a vll )/L li L v , (1) 

where L M , /i = 1 , . . . , 4 are the linear extensions of 
the hypertorus, 

<V, - a Vfl = /Nl - B W /Nk . (2) 
The integers n$ and n$ summed to n^ v — n^} + 

(2) 

riftJ define the so-called twist tensor n^ v . For n^ v = 
(mod. N) the twists are commuting and can be 
removed by appropriate gauge transformations such 
that gauge fields become periodic. 

For n X 2 = —n^2 = 1 (with other components equal 
to zero) n^„ = 0, oi\2 — &2\ = 1/kl we get a constant 
magnetic field in the third direction B 3 = Fyi- The 
action of this field on the hypertorus with L\ = L2 = 
L3 = L s and L4 = L t is equal to 

S DS = l/2g 2 (B$) 2 V 4 = l/g 2 Tr(B 3 ) 2 V 4 

= 87r 2 /ff 2 x N/2M x L t /L s . (3) 

Thus, for 517(2) S DS = S mst L t /L s , for 517(3) S DS = 
S/ASinstLt/Ls, where the instanton action is Si ns t = 
8n 2 /g 2 . In the SU{2) case the magnetic field B 3 is 
equal to B 3 = 2ir diag(l, — 1)/L 2 and its flux $ over 
the 12-plane of the hypertorus is a multiple of 27r: 
$ = 2ir diag(l, — 1) . This means that in the periodic 
gauge such a field could remain Abelian because of 
exp(i<£>) = 1 . In the SU (3) case the magnetic field 
B 3 = 7r diag(l, 1, —2)/L 2 has a flux over the 12-plane 
of the hypertorus equal to $ = 7r diag(l, 1,-2). Now 
exp(i$) = diag(— 1, — 1, 1) is not equal to the unity 
matrix and this means that in the periodic gauge such 
a field could not remain Abelian. 



III. 577(2) EMBEDDED DIRAC SHEET 
SOLUTIONS 

The Dirac sheet seen on the lattice in the SU(2) 
case [TH, [lH is observed also in SU(3) lattice sim- 
ulations. We will call it standard DS. New, specific 
for the SU(3) case, is the Dirac sheet with an action 
value equal to 3/4 of the action of the standard DS. 
In the following we will call it non-standard DS. It is 
also seen in lattice simulations. 

In SU(2) a constant Abelian magnetic field is not 
stable under fluctuations of the gauge field. Charged 
(off diagonal) components of the gauge field have a 
Savvidy eigenmode [19( with negative eigenvalue 

A = -4n/L 2 s . (4) 



The situation can be stabilized by introducing a con- 
stant Abelian scalar potential A\. Normally a con- 
stant Abelian scalar potential can be gauged away. In 
our case due to periodicity in time direction it can be 
gauged away only modulo 2it/L t . The interaction of 
charged (off-diagonal) components of the gauge field 
with this potential adds a positive term to the eigen- 
value A, turning it into 

A = -4n/L 2 s + (A\f . (5) 

The presence of the scalar potential leads to a non- 
trivial holonomy H that is defined as 

H= lim P cxp(i / A 4 (S,t)dt). (6) 

The holonomy is parametrized as H = 
diag(e 27r ^ 1 ,e 2 ' riA12 ) with m < fi 2 < M3 = 1 + Mi and 
Ml + M2 = 0. Thus, positive numbers mi = ^2 — Mi> 
m 2 = M3 — M2 sum up to unity mi + TO2 = 1. The 
eigenvalue A then becomes equal to 

A = -4tt/L 2 + (2 7 rmi/ J L t ) 2 , (7) 

and its positiveness requires L t /L s ^/Tr < mi 2 < 
1 — L t j L s ^pK . Therefore, nontrivial holonomy stabi- 
lizes DS and just this situation was observed in SU(2) 
lattice cooling .12] and elucidated in Ref. [Hj]. 

Now let us consider the embedding of this standard 
DS event into SU(3) group. Let vector potentials A\ 2 
be proportional to diag(l, — 1, 0) and the scalar poten- 
tial to give the holonomy 

H = diag(e 2rrvi , e 27 ^ 2 , e 27 "" 3 ) (8) 

with Hi < \i 2 < M3 < tM = 1+Mi an d M1+M2+M3 = 0. 
Now three positive numbers mi = /J2 — Mi) m 2 = 
M3— M2 ! m 3 = M4 — M3 sum to unity mi +m2 + m3 = 1. 
Stability of the DS under fluctuations of charged (off- 
diagonal) (1,2) — (2, 1) components of the gauge fields 
requires Lt/L s y / 7r < mi < 1 — L t /L s y / 7r. The other 
off-diagonal (2,3) — (3,2) and (3,1) — (1,3) compo- 
nents of the gauge fields have charges with respect to 
the diag(l,— 1,0) generator of the SU(3) group be- 
ing two times smaller than the (1,2) — (2,1) com- 
ponents. Hence the stability of DS under their fluc- 
tuations requires LtjL s \f2/K < mi < 1 — Lt/L s ^/2w 
and Lt/L s y/2~TT < m 3 < 1 — Lt/ L s ^/2~k, correspond- 
ingly. Taking into account that the magnetic Abelian 
field could lie also in other SU{2) subgroups of the 
5f7(3) group, i.e. would then be proportional to 
diag(l,0, — 1) or to diag(0,l,— 1) generators, we see 
that the standard DS in SU(3) group will be stable 
for values of the holonomy restricted by the following 
constraints on the holonomy parameters 7711,7712,7713 

L t /L s V2^<m h2 . 3 <l-L t /L s V2^. (9) 

We shall visualize the stability criteria in a (X, Y) 
plot in the complex plane, X — 5f(l/3 Tr H) and 
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FIG. 1: The S?7 (3) triangle and the inscribed region of sta- 
bility expected for standard DS configurations (enclosed 
by the dashed line) compared with standard DS events 
found in actual lattice cooling (filled circles). 



Y = 9(1/3 Tr H)). The corresponding region for the 
standard DS configurations is shown on Fig. 1. The 
external curved triangle encloses all possible values 
of one third of the trace of an unitary matrix (the 
holonomy) that can be obtained by the variation of 
the phase parameters mi,m,2,m3 in the region < 
"ti,2,3 < 1, while the sum is constrained by mi + 
m, 2 + m 3 = 1. The smaller, inscribed curved triangle 
(bounded by the dashed line) is the region of stability 
of standard DS events. 



IV. NON-STANDARD DIRAC SHEETS 

Coming now to the discussion of the stability of 
non-standard DS solutions one should first stress 
that by construction constant Abelian magnetic fields 
can be supplemented only by a constant Abelian 
scalar potential proportional to the same diagonal 
SU(3) generator to which the magnetic field is pro- 
portional. If the magnetic field is equal to B3 = 
7T diag(l, 1, —2)/Lg, then in a constant Abelian scalar 
potential 

At = diag(27r^ 1 /L t , 2nn 2 /L t , 2vr/i 3 /£ t ) (10) 

the holonomy parameters [i\ and \ii should be equal 
to each other: \i\ = /12 (mi = ). The fluctuations 
of the (1, 2) — (2, 1) components of gauge fields in this 
case do not interact with both the magnetic field and 
the static scalar potential. For fluctuations of charged 
(2, 3) — (3, 2) and (3, 1) — (1, 3) components the lowest 
modes have eigenvalues 

A 23 - -3n/L 2 s + (2irm 2 /L t ) 2 (11) 

and 

A13 = -3n/Ll + (2irm 3 /L t ) 2 (12) 
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FIG. 2: The region of stability of non-standard DS config- 
urations (the three sides of the unclosed SU(3) triangle) 
compared with non-standard DS events found in actual 
lattice cooling (filled circles). 

correspondingly. So, the stability of such non- 
standard DS solutions is possible for 

mi = 0, y/3/4nL t /L s < m 2 , 3 < 1 - y^/4nL t /L s . 

(13) 

For other non-standard DS solutions the region 
of stability can be obtained by the permutations of 
holonomy parameters mi, m2, ma- The stability re- 
gion is shown in the (A, Y) plot of Fig. 2 and happens 
to coincide with the boundary of the unclosed SU(3) 
triangle of Fig. 1. 

V. NUMERICAL RESULTS 

For a numerical study of standard and non-standard 
DS solutions we have employed the standard Wil- 
son plaquette action Sw, creating an ensemble with 
/3 = 6/g 2 where g denotes the bare coupling constant. 
On a lattice for Lt = 4, L s = 16 the coupling constant 
related to the first order deconfinement transition is 
equal to j3d ~ 5.69. The initial Monte Carlo ensem- 
ble was generated in the confined phase at (3 = 5.63. 
As expected, this has guaranteed that in the process 
of cooling the holonomy has remained sufficiently non- 
trivial, such that the emerging DS configurations were 
stable. We have found configurations stable against 
further cooling with the action S = 1/ASinst and 
S = 3/\QSinst in perfect agreement with analytical 
knowledge. We have stopped cooling at the moment, 
when the relative variation of action density inside 
the configuration became smaller than 10~ 4 (homoge- 
neous configurations) and have measured the value of 
holonomy (the average Polyakov loop) . The Polyakov 
loop also has happened homogeneous. The distance 
of local values of it from the average value was not 
larger than 10~ 5 . The scatter plots of DS events in 
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the (X, Y) plane of the real and imaginary part of 
the Polyakov loop are shown in Figs. 1 and 2. The 
dots lie perfectly inside the regions of stability for the 
respective type of DS configurations. The configura- 
tions obtained turned out to be purely magnetic and 
- applying maximally Abelian gauge - show constant 
Abelian magnetic fluxes. 

We did not particularly attempt to find Dirac sheets 
at higher temperature, f3 > f3 c . We know from other 
simulations that the holonomy of such equilibrium 
configurations under cooling rapidly evolves towards 
central elements where Dirac sheets are unstable and 
therefore would have escaped observation. 

VI. CONCLUSION 

In conclusion, purely Abelian constant magnetic 
field configurations have been observed emerging from 
the process of cooling equilibrium (Monte Carlo) lat- 



tice fields representing the confined phase of SU(3) 
gluodynamics. They were found to be absolutely sta- 
ble provided their Polyakov loop was sufficiently non- 
trivial. We have shown here that this fact is related 
to the notion of marginal stability of the appropriate 
constant magnetic field configurations. 

Finally we have to admit that the Dirac sheet con- 
figurations discussed in this paper will not play any 
role in the thermodynamic limit of the theory since 
their action tends to zero in this limit. 
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